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The dominant X-ray wind in massive star binaries 
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Abstract. We investigate which shocked wind is responsible for the majority of the X-ray emission in colliding 
wind binaries, an issue where there is some confusion in the literature, and which we show is more complicated 
than has been assumed. We find that where both winds rapidly cool (typically close binaries) , the ratio of the wind 
speeds is often more important than the momentum ratio, because it controls the energy flux ratio, and the faster 
wind is generally the dominant emitter. When both winds are largely adiabatic (typically long-period binaries), 
the slower and denser wind will cool faster and the stronger wind generally dominates the X-ray luminosity. 
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1. Introduction 

The violent wind-wind collision in massive star binaries 
creates a region of high temperature shock-heated plasma, 
which can contribute to the total system emission at ra- 
dio, infrared, optical, ultraviolet and X-ray wavelengths. 
Over the last 20 years, theoretical models have focussed 
mainly on the dynamics of the stellar winds and the wind 
collision zone (WCZ), and on the resulting X-ray emis- 
sion (Pittard ^000| and references therein). As the X-ray 
emission is dependent on the physical conditions within 
the WCZ and on the distribution and properties of the 
unshocked attenuating wind material, observations pro- 
vide information on basic parameters of the system (see 
Stevens et al.|1996|: Zhekov & Skinner 120001 Pittard & 



Corcoran 2002) 



To interpret such observations, it is useful to know 
which wind dominates the X-ray emission. For instance, if 
there is a velocity difference between the winds, we would 
like to be able to predict how hot the broad-band emission 
will be, or if there are abundance differences, how strong 
or weak certain lines will be (this is of particular relevance 
now that grating observations can resolve line profiles). 

Simple expressions for the X-ray luminosity from each 
of the shocked winds were presented in a complex, yet el- 
egant paper, which provides some answers to these ques- 
tions in terms of fundamental parameters of the system 
(Usov 1992). Their applicability has resulted in their com- 
mon use in the literature and in observing proposals. 
However, we have recently discovered inconsistencies be- 



tween results from these expressions and those determined 
from numerical models, leading to some confusion on the 
issue of the dominant X-ray emitting wind. In this letter 
we reinvestigate previous conclusions in the literature. 

2. Estimates of L1/L2 

Let us define Li as the X-ray luminosity from the shocked 
wind with the greater momentum flux (i.e. the 'stronger' 
wind), and L2 the equivalent from the shocked weaker 
wind. Analytical estimates in the literature of the ratio 
of L1/L2 exist for two limiting cases: the radiative limit, 
where the cooling timescale for the hot gas from both 
winds is assumed small in comparison to the timescale for 
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flow of this gas out of the system (i.e. x < Ij see Sec. 2.2); 
and the adiabatic limit, where the opposite is true. We 
will first focus on these limiting cases, before discussing 
the behaviour of L1/L2 between these limits. 

2.1. L1/L2 in the radiative limit 

In this limit, the entire kinetic energy thermalized by the 
shocks is immediately radiated (normally with the major- 
ity at X-ray energies), and the region of shocked gas is 
thin. Binary systems with a small orbital separation are 
favoured, and examples near this limit include V444 Cyg 
(WN5 + 06;P = 4.2 d; Corcoran et al.|l99§) a nd DH Cep 
(05.5-6V + 06.5-7V; P = 2.11 d; Penny et aL p^Z] ). 

Let us first consider the region of the WCZ which lies 
directly between the stars. As there is a momentum bal- 
ance, the relative kinetic energy fluxes are proportional to 
the ratio of the wind speeds, so emission from this vol- 
ume should be dominated by gas from the star with the 



2 



Pittard & Stevens: The dominant X-ray wind in massive star binaries 



faster wind. If, for simplicity, we initially assume that the 
pre-shock wind speeds are spatially invariant, we find that 
the majority of the thermalized energy in the WCZ occurs 
close to the line of centres (see Fig. ^ . Under these con- 
ditions, the faster wind dominates the total emission. 

If the pre-shock wind speeds are approximately equal, 
neither star should significantly dominate the total emis- 
sion. In such circumstances, one might expect secondary 
considerations, such as the wind momentum ratio, to be- 
come important. Luo et al. ( 199C| ) noted that if vi ~ V2, 
the weaker wind should be more luminous, based on the 
premise that it has, on average, a larger velocity compo- 
nent normal to the shock front: referring to V444 Cyg 
they argue that, "Because of this, the shocked 06 wind 
dominates the X-ray emission from the shocked region" . 
However, we show below that even for very large wind 
momentum ratios, the stronger wind still accounts for al- 
most half of the total emission. Therefore Luo et al. 's use 
of the word "dominates" goes too far, although their ba- 
sic conclusion that the 06 wind is the majority emitter 
is consistent with our analysis. Complicating factors, such 
as wind acceleration, are discussed later. 

We can obtain analytical es timat es for L1/L2 from 
the relevant equations in Usov ( |l992|) . For the radiative 
limit, the total X-ray luminosity from the external wind 
shock is (cf. his Eq. 88) Li = 2.3 poo vf Zqq, where 
poo = Mi/4:TtD'^vi, zqb is the distance of the contact 
discontinuity from the centre of the star with the weaker 
wind (= Y^D/(1 -I- y^)), D is the stellar separation, and 
rj is t he mo mentum ratio of the winds (= M2V2/ MiVi). 
Usov ( |l992[) parameterizes the luminosity from the inter- 
nal wind shock as L2 = 0.22M2f|. 

With the same parameter values as used by Luo et 
al. ( |1990| ) for V444 Cyg, these equations yield L1/L2 = 
0.59 (i.e. that the weaker wind is the dominant X-ray emit- 
ter), in agreement with the statement in Luo et al. (19901). 



From Usov (1992) we find that L1/L2 oc Miwf /Afaul 



pivf/ P2V2, which at the stagnation point is proportional 
to Vi/v2, such that as expected the faster wind will nor- 
mally be the dominant X-ray emitter. While this agree- 
ment is satisfying, an i ndicat ion of potential problems with 
the equations in Usov ( 1992 ) is revealed by considering the 
predicted value of L1/L2 for identical winds: in this case 
we obtain L1/L2 = 0.2, instead of unity! 

We have therefore performed a numerical calculation 
of the kinetic energy flux normal to the wind shocks to 
precisely determine L1/L2 for various wind parameters. 
The position of the contact discontinuity (which is also the 
position of the wind shocks since the wind collision region 
is thin) was obtained from integration of Eq. 4 in Stevens 
et al. (1992). We find that the value of L1/L2 primarily 
depends on the ratio of the wind speeds and to a much 
lesser extent on the wind momentum ratio. Results are 
shown in Fig. ^ and in Table ^ where we have defined 
S to be the fractional wind kinetic power normal to the 
contact discontinuity i.e. Li ~ O.SSiAfiu^ for the stronger 
wind. 



For 77 = (0.01, 0.1, 1.0), we find that Si = 
(0.0042, 0.033, 0.167) and S2 = (0.564, 0.403, 0.167) 
respectively (the latter values being the analytical limit of 
1/6), consistent with the solid angle of the wind collision 
zone as viewed from the star with the stronger (weaker) 
wind decreasing (increasing) with decreasing 77. The ra- 
tio of Usov's equations yield values for L1/L2 which are 
too low in comparison to the exact numerical calculation 
by factors of (1.09, 1.74, 4.88) for rj = (0.01, 0.1, 1) re- 
spectively, irrespective of the ratio of vi /v2 ■ Thus Usov's 
equations for the radiative limit are somewhat in error, 
being most accurate for low values of 77. 

Fig. 1^ also shows the gradient of S as a function of 
the off-axis distance, r. In Table ^ we list the value of r 
at which J'J^g(dS/dr)dr is 50% and 90% of the asymp- 
totic value of S. For equal strength winds, 50% and 90% 
of the maximum wind kinetic power thermalized by the 
shocks (and hence radiated) occurs within r/D = 0.38 and 
0.95 (this off-axis distance is much smaller than the cor- 
responding value for adiabatic winds cf. Luo et al. 1990). 
Therefore distortion of the WCZ by orbital motion should 
not significantly alter our previous conclusions. We can 
also safely relax our assumption of spatially invariant wind 
speeds as we only need to consider this ratio between the 
stars. 

Finally, for close binaries we should also consider the 
potential role of radiative acceleration, inhibition (Stevens 
& Pollock |1994D and braking (Gayley et al. |1997D of the 
winds on L1/L2. As already shown, the ratio of the wind 
speeds is the critical parameter in such systems, and it 
seems plausible to expect this to continue to be true in 
the presence of these 3 effects. In such binaries, it is prob- 
able that neither wind will have room to reach terminal 
speed. Since the stronger wind will have more room to 
accelerate, we would normally expect it to be faster, and 
to dominate the emission. However, if there is substan- 
tial radiative inhibition plus braking, this will cause the 
strong wind to be by far the slower and will likely shift 
the dominant energy generation back to the weaker wind. 

2.2. L1/L2 in the adiabatic limit 

In direct contrast to the radiative limit, colliding wind 
systems are most likely to be near the adiabatic limit if 
the stellar separation is large. Perhaps the best example of 
such a system is the Wolf-Rayet binary WR 140 (WC7 -I- 
04V; P = 2900 d; WiUiams et al. [19971) . Luo et al. ( |1990D 
argue that for a WR-l-0 system in the adiabatic limit "the 
shocked WR stellar wind dominates the X-ray emission." 



Myasnikov & Zhekov (1993) arrive at a similar conclusion 
for their "standard system" which is also near this limit: 
"the whole of the luminosity is due entirely to the emission 
of the shocked gas of the WR wind" , being more than an 
order of magnitude greater than the emission from the 
weaker wind.^ 



^ Although Myasnikov & Zhekov ( 199S ) actually use different 
abundances for each wind, gas with typical WN abundances 
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Table 1. Ratio of the wind X-ray luminosities in the ra- 
diative hmit assuming spatially invariant pre-shock wind 
speeds. The value of Li/ L2 primarily depends on the ratio 
of the wind speeds, Wi/w2, and to a lesser degree on the 
wind momentum ratio, 77. Also tabulated are the ratios of 
the mass-loss rates and the wind kinetic power. 
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Table 2. Values of the off-axis distance, r, at which Si and 
S2 are 50% and 90% of their maximum value for various 
77, and assuming spatially invariant, equal speed winds. 
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Fig. 1. Results for a wind- wind collision in the radiative 
limit, with equal and spatially invariant pre-shock speeds. 
The top panel shows the position of the contact disconti- 
nuity as a function of wind momentum ratio, 77. Star 1 is 
located at (r,z) = (0,0), and star 2 at (r,z) = (0,1). The 
middle and bottom panels show the value of dS/dr (see 
text) as a function of r and r/ for the wind of star 1 and 
star 2 respectively. Note that the majority of the thermal- 
ized kinetic power of the winds occurs well within r — D 
(see also Table ||). 



While the appropriate equations in Usov (1992) are 
known to be lower limits (due to the omission of line emis- 
sion), the ratio of L1/L2 should not be overly affected and 
can once again be calculated. Using the same values for the 
stellar wind parameters as Myasnikov & Zhekov ( 1993D , 
Eqs. 89 and 95 in Usov ( |1992D yield L1/L2 = 1-36, indi- 
cating that neither wind is particularly dominant. This is 
clearly in disagr eement with the published statem ents in 
Luo et al. ( |l990|) and Myasnikov & Zhekov (|l993|) . 



has a similar emissivity as gas with solar abundances, such 
that this difference is unimportant to their findings. 



To investigate this issue we have calculated the X- 
ray luminosity from numerical simulations of the WCZ 
near the adiabatic limit. Table ^ summarizes our findings 
for winds with equal, spatially invariant pre-shock speeds. 
The general trend is for L1/L2 to increase with decreasing 
77 (which is opposite to the radiative limit), and for this 
ratio to become very large for sm all 77. W e thus find agree- 
ment wit h the work of Luo et al. (1990) and Myasnikov & 
Zhekov (1993), and conclude that the appropriate equa- 
tions in Usov ( 1992 ) are again somewhat in error. 

The underlying reason for the trend shown in Table || 
concerns the ratio of the cooling timescale to the flow 
timescale in each of the winds. Stevens et al. ( 1992| ) noted 
that, near the local minimum in the cooling curve (which 
implies that the pre-shock velocity at the stagnation point 
is in the range 900kms~^ ^ v < 3500 kms^^), this ra- 
tio can be approximated as x ~ '^'l di2/M^T, where vs is 
the wind velocity in units of 1000 km s~^ , di2 is the dis- 
tance from the star to the contact discontinuity in units 
of 10^^ cm, and M_7 is the mass-loss rate of the star 
in units of lO~^M0yr~-'^ . The shocked wind is radiative 
when X < 1, and approaches the adiabatic limit when 
X S> 1. If we ignore the dependence of x on d (i.e. if we 
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Table 3. Ratio of the wind X-ray luminosities as a func- 
tion of wind momentum ratio, 77, for colliding wind sys- 
tems with equal, spatially invariant wind speeds, equal 
abundances, and near to the adiabatic limit. 
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were to assume that the relevant distance appropriate to 
the flow dynamics is the same for each wind), the ratio of 
this characteristic cooling parameter for the two winds is 
X1/X2 ^ M2vt/Mivl 

For the results in Table || where vi = V2, X1/X2 ~ 
M2/M1 = r). Hence as the value of r] decreases in Table ||, 
the stronger wind becomes more efficient at emitting X- 
rays relative to the weaker wind. This is consistent with 
the fact that although the post-shock density at the stag- 
nation point is the same for both winds, its decline with 
off-axis distance is faster for th e wea ker wind (as observed 
also by Myasnikov & Zhekov 1993| ). We note also that 
within a given distance from the stagnation point, the vol- 
ume occupied by the shocked stronger wind exceeds that 
occupied by the shocked weaker wind. 

One might expect that L1/L2 1 when Xi/X2 = 1, ir- 
respective of the value of 77. Again ignoring the dependence 
of X on d, we find that rj = {v2/vi)^ for X1/X2 = 1- Thus 
to obtain comparable luminosity from each shocked wind 
when 77 = 0.1, we require vi ^ 1.6 V2- Since, in reality, 
mass-loss rates from early-type stars can vary by several 
orders of magnitude, whereas wind velocities lie typically 
within the range 1000- 3000 kms'^ (excluding LBV's), 
we expect most systems near the adiabatic limit will have 
a luminosity dominated by the stronger wind. 

2.3. L1/L2 in-between these limits 

Estimates of the luminosity ratio for systems where the 
shocked winds are in-between the limiting radiative and 
adiabatic cases must be done on a case-by-case basis. 
Somewhere in this region of parameter space, the wind 
which dominates the X-ray emission must switch over from 
the one with the faster pre-shock speed (radiative limit) to 
the one with the slower pre-shock speed (adiabatic limit). 
In contrast, situations where one wind is clearly radiative 
and the other is closer to being adiabatic will have their X- 
ray emission dominated by the former (e.g. in 7^ Velorum, 
where xi 1 s-nd X2 > 1, the X-ray emission is domi- 
nated by the shocked WR wind) . 

3. Summary 

In this letter we have re-examined the issue of which wind 
is the dominant X-ray emitter in colliding wind binaries. 



following the discovery of some confusion in the existing 
literature. Our work supports the earlier conclusions in 
Luo et al. ( |199GD and Myasnikov & Zhekov ( |1993D , though 
is sometimes in disagreement with the relevant analytical 
equations in Usov (1992). 

For systems near the radiative limit (typically short pe- 
riod binaries), we find that the primary infiuence on L1/L2 
is the ratio of the wind speeds, vi/v2, since it controls the 
energy flux ratio. For vi — V2t Li ^ L2 irrespective of the 
wind momentum ratio, 77. When vi 7^ U2, the faster wind 
normally dominates the luminosity, although it is unlikely 
t o do so by more than a factor of 5. The equations in Usov 
(1992) predict values for L1/L2 which can be in error by 
up to a factor of 4. 

For systems near the adiabatic limit (e.g. long-period, 
high eccentricity binaries at apastron), we confirm earlier 
findings that the stronger wind is typically the dominant 
X-ray emitter, often by an order of magnitude relative to 
the weaker wind. This is because the dominant driver of 
L1/L2 is the ratio of Xi/X2: so that the stronger wind 
can in general more easily radiatively cool. If we force 
Xi/X2 ~ Ij we find that both shocked winds contribute 
roughly equally to the X-ray emission irrespective of the 
value of 77. In contrast, the equations in Usov ( |1992| ) yield 
Li ~ L2 irrespective of the assumed wind parameters. 

For systems in-between these limits, we anticipate that 
the dominant luminosity should generally switch from the 
faster wind (radiative limit) to the slower wind (adiabatic 
limit). Where one shocked wind is substantially closer to 
the radiative limit than that of the other, the X-ray emis- 
sion will naturally be dominated by the former. 

This work provides a basic understanding of the dom- 
inant factors controlling the luminosity ratio in colliding 
wind binaries. Detailed calculations are needed to further 
investigate the effects of wind acceleration/braking and 
different wind abundances. 
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